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The properties of two-flavored massive Nambu-Jona-Lasinio model in (l+l)-dimensional space- 
time are investigated in the presence of isospin and quark number chemical potentials. The consid- 
eration is performed in the large- A'^c limit, where Nc is the number of colored quarks. It is shown in 
the framework of this model that charged pion condensation phenomenon of dense quark/hadron 
isotopically asymmetric matter is rather a spatially inhomogeneous than a homogeneous one. 



I. INTRODUCTION 
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^jr^, Recently, much attention has been paid to the investigation of the QCD phase diagram in terms of baryonic as well 
^ ■ as isotopic (isospin) chemical potentials. The reason is that dense baryonic matter which can appear in heavy-ion 
collision experiments has an evident isospin asymmetry. Moreover, the dense hadronic/quark matter inside compact 
stars is also expected to be isotopically asymmetric. To describe the above mentioned realistic situations, i.e. when 
the baryonic density is comparatively low, usually different nonperturbative methods or effective theories such as 
chiral effective Lagrangians and especially Nambu - Jona-Lasinio (NJL) type models [l| are employed. In this way, 
I— 1' the QCD phase diagram including chiral symmetry restoration color superconductivity 0t11j ^^.^ charged 

', pion condensation (PC) phenomena [Tol - [l^ were investigated under heavy- ion experimental and/or compact star 
conditions, i.e. in the presence of such external conditions as temperature, chemical potentials and possible external 
Q_^^. (chromo)magnetic fields (see the above references). 

(]J • Obviously, the (3-1-1) -dimensional NJL models depend on the cutoff parameter which is typically chosen to be of 
' the order of 1 GeV, so that the results of their usage are valid only at comparatively low energies, temperatures and 
densities (chemical potentials). Moreover, there exists also a class of renormalizable theories, the (l-l-l)-dimensional 
' chiral Gross-Neveu (GN) type models [l^ [13 ; that can be used as a laboratory for the qualitative simulation of 
^ , specific properties of QCD at arbitrary energies. Renormalizability, asymptotic freedom, as well as the spontaneous 
^\ breaking of chiral symmetry (in vacuum) are the most fundamental inherent features both for QCD and all GN type 
models. In addition, the /i — T phase diagram is qualitatively the same for the QCD and GN models [l7l - [20l | (here /x is 
the quark number chemical potential and T is the temperature). Note also that the GN type models are suitable for 
the description of physics in quasi one-dimensional condensed matter systems like polyacetylene (2]| . It is currently 
well understood (see, e.g., the discussion in [l^, H^, [13) that the usual no-go theorem ^2^, which generally forbids the 
spontaneous breaking of any continuous symmetry in two-dimensional spacetime does not work in the limit Nc — > 00, 
, where Nc is the number of colored quarks. This follows from the fact that in the limit of large Nc the quantum 
. . ■ fluctuations, which would otherwise destroy a long-range order corresponding to a spontaneous symmetry breaking, 
^ [ are suppressed by 1/Nc factors. Thus, the effects inherent for real dense quark matter, such as Cooper pairing 
phenomenon (spontaneous breaking of the continuous U{1) symmetry) or charged pion condensation (spontaneous 
breaking of the continuous isospin symmetry) might be simulated in terms of a simpler (l+l)-dimcnsional GN-type 
, model, though only in the leading order of the large Nc approximation (see, e.g., 20, 24] and [25 27, 30], respectively). 
This paper is devoted to investigation of the charged pion condensation (PC) phenomenon in the framework of 
the (l-l-l)-dimensional NJL model with two quark flavors and in the presence of the quark number (p) as well as 
isospin (/X/) chemical potentials. The consideration is performed in the leading order of the 1/A^c-expansion. In our 
previous papers |25l - l27| the phase diagram of the above mentioned massless or massive NJL2 model was already 
investigated in the case of homogeneous, i.e. independent of space coordinate, order parameters (chiral and charged 
pion condensates). The situation corresponds to the conserved Lorentz and spatial translational invariance and is 
adequate to physical systems in vacuum, i.e. at zero chemical potentials. However, in dense baryonic matter, i.e. 
at nonzero quark number chemical potential, there might appear new phases with a spatially inhomogeneous chiral 
and/or charged pion condensates which destroy both chiral and/or isospin as well as spatial translational invariances 
of a system. In particular, the possibility of the phase with inhomogeneous chiral condensate was discussed in the 
framework of both (l-l-l)-dimensional [22, 28 30] and (3-|-l)-dimensional [3ll - l38| models. At the same time the 
possibility of the phase with inhomogeneous charged pion condensation is less investigated. 

Thus, in this paper, in contrast to [25l - [27j . we consider the phase portrait of the above mentioned massive (1+1)- 
dimensional NJL model with two chemical potentials, fi and /x/, in the leading order of the 1/A^c-expansion taking into 



^ Below we shall use the notation "NJL2 model" instead of "chiral GN model" for (H-l)-dimensional models with a continuous chiral 
symmetry, since the chiral structure of the Lagrangian is the same as that of the (3+l)-dimensional NJL model. 

^ For example, spatially inhomogeneous ansatz for the charged pion condensate was investigated in the framework of the two-flavored 
NJL4 model in [T3| . It was claimed there that inhomogeneous PC phase is possible only at rather high values of an isotopic chemical 
potential, /ii > A, where A ~ 0.65 MeV is a cutoff parameter. So this result might be out of the scope of a model application. Moreover, 
the authors of I13l| made some technical simplifications in their research of the inhomogeneous PC phenomenon. 
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account the possibility that the charged pion condensate might become spatially inhomogeneous. The temperature 
is taken to be zero. For simplicity, for the chiral condensate we use a spatially homogeneous ansatz. (In contrast, in 
our previous paper [30l | the spatially inhomogeneous ansatz for the chiral condensate and homogeneous one for the 
charged pion condensate was used in the framework of the same massless NJL2 model.) Notice once more, the isotopic 
asymmetry is an inevitable property of dense quark matter which might be created in heavy-ion collision experiments 
or inside compact stars. So, we believe that such a simplified study in the framework of the two-dimensional NJL 
model with isospin chemical potential could shed new light on the properties of real dense baryonic matter and will 
provide a deeper understanding of the charged PC phenomenon. In particular, it is shown in our paper that in the 
framework of the model under consideration a PC phase with nonzero baryon density is realized just in the case 
of inhomogeneous charged pion condensate but not in the case of homogeneous one. In analogy, one can expect 
that in real (3-|-l)-dimensional dense hadronic/quark matter the charged PC phenomenon is realized rather with 
inhomogeneous pion condensate than with spatially homogeneous one. 

The paper is organized as follows. In Section II we derive, in the leading order of the large iVc-expansion, the 
general expression for the thermodynamic potential of the two-flavored massive NJL2 model with quark number 
chemical potential /i and isospin chemical potential fij in the case of spatially homogeneous chiral condensate and 
inhomogeneous PC. First, in Sec. Ill we reduce our consideration to the case of homogeneous PC and find that in this 
case only a PC phase with zero density of quarks is possible in the model. Second, in Sec. IV it is shown that charged 
PC phase with nonzero quark density in the framework of the model is possible only with spatially inhomogeneous 
charged pion condensate. Final Sec. V presents a summary and some concluding remarks. The discussion of some 
technical problems are relegated to two Appendices. 



II. THE MODEL AND ITS THERMODYNAMIC POTENTIAL 



We consider a (l-l-l)-dimensional NJL2 model to mimic the phase structure of real dense quark matter composed 
of two massive quark flavors (u- and d- quarks) . Its Lagrangian has the form: 



L 



Yidp - mo + fij" + yT37*^ 



G 



{qqf + [qh^fqf 



(1) 



where the quark field q{x) = qia{x) is a flavor doublet (i = 1, 2 or i = m, d) and color A^c-plet (a = 1, Nc) as well as 
a two-component Dirac spinor (the summation in ([T]) over flavor, color, and spinor indices is implied); (fc = 1, 2,3) 
are Pauli matrices; the quark number chemical potential in ([1]) is responsible for the nonzero baryonic density of 
quark matter, whereas the isospin chemical potential /i/ is taken into account in order to study properties of quark 
matter at nonzero isospin densities (in this case the densities of u and d quarks are different). The Dirac gamma 
matrices in two-dimensional spacetime have the following form: 
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(2) 



Evidently, the model ([T|) is a simple generalization of the original (l-l-l)-dimensional Gross-Neveu model [Ts!] with 
a single massless quark color A^c-plet to the case of two massive quark flavors and additional chemical potentials. 
As a result, in the case under consideration we have a modifled flavor symmetry group, which depends essentially 
on whether the bare quark mass ttiq and isospin chemical potential /i/ take zero or nonzero values. Indeed, at 
fii = 0,TOo — the Lagrangian ([T]) is invariant under transformations from the chiral SUl{'2) x SUii{2) group. Then, 
at /i/ 7^ 0,TOo = this symmetry is reduced to Ui^l{^) x Ui^ji{l), where = T3/2 is the third component of the 
isospin operator (here and above the subscripts L, R mean that the corresponding group acts only on left, right handed 
spinors, respectively). Evidently, this symmetry can also be presented as Ui^{l) x [7^/3(1), where Ui.^{l) is the isospin 
subgroup and Uai^ (1) is the axial isospin subgroup. Quarks are transformed under these subgroups as g — exp{iaT3)q 
and q — exp(ia7^r3)g, respectively. In the case toq 7^ 0,/j,/ = the Lagrangian ([T]) is invariant with respect to the 
SUi{2), which is a diagonal subgroup of the chiral SUl{2) x SUr{2) group. Finally, in the most general case with 
Too 7^ 0,/L(/ 7^ the initial model ([1]) is symmetric under the above mentioned isospin subgroup Ui^{l). In addition, 
in all foregoing cases the model is invariant under color SU(A^c)-, baryon charge J7_b(1)- and electric charge Uq{1) 
groups. 

The linearized version of the Lagrangian ([1]), which contains composite bosonic flclds a{x) and "Kaix) (a = 1, 2, 3), 
has the following form: 



Y'^dp - Too -h M7° + ■^'^37° - Cr - l-f^lTaTa 



El 

4G 



aa + n„TT„ 



From the Lagrangian ([3]) one obtains the following constraint equations for the bosonic flelds 

G G 
f^ix) = -2 — {qq); TTaix) = -2 — {qij^Taq). 



(3) 



(4) 
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Obviously, the Lagrangian ([3]) is equivalent to the Lagrangian ([T]) when using the constraint equations (|4]). Fur- 
thermore, it is clear that the bosonic fields ^ are transformed under the isospin Ui^{l) subgroup in the following 
manner: 

[7/3(1) : (7 — > cr; TTa — > tt^; tti — > cos(2a)7ri + sin(2a)7r2; 7r2 cos(2a)7r2 — sin(2a)7ri, (5) 

i.e the expression (ttJ + 7r|) remains unchanged under transformations of the isospin subgroup Uj^{l). 

To avoid the no-go theorem (23j . which forbids the spontaneous breaking of continuous symmetries in the considered 
case of one space dimension, we restrict the discussion only to the leading order of the large Nc expansion (i.e. to the 
case Nc — 00), where this theorem is not valid [l^, [10, HI] • In particular, the effective action Ses{<^,T^a) can be found 
in this approximation through the relation: 



exp(i5cff(cr,7ra)) = N' J [dq][dq] exp (^i J Ld'^x^, (6) 
where N' is a normalization constant. It is clear from ([3|) and ^ that 



5cff (a, TTa) = -Nc J "- d'x + 5eff , (7) 

where the quark contribution to the effective action, i.e. the term ScS in (O, is given as follows 

exp(i5efF) = N' J[dq][dq] exp(i J q[iYdp - toq + M7° + t^Tsl" - cr - ij^TTaTalq d^x^ . (8) 

Here we used the notation u = fii/2. The ground state expectation values (cr(x)) and (iTaix)) of the composite bosonic 
fields are determined by the saddle point equations. 



Sa{x) ' <^7ra(a;) 

where a = 1,2,3. In vacuum, i.e. in the state corresponding to an empty space with zero particle density and zero 
values of the chemical potentials fi and /i/, the quantities (cr(x)) and (7ra(x)) do not depend on space coordinates. 
However, in dense quark medium, when /i 7^ 0, /i/ 7^ 0, the ground state expectation values of bosonic fields might 
have a nontrivial dependence on x. In particular, in this paper we will use the following ansatz: 

{a{x))=M-ma, (773(2;)) = 0, (7ri(x)) = Acos(25a;), (772(0;)) = A sin(25a;), (10) 

where M, b, and A are constant quantities. In fact, they are coordinates of the global minimum point of the thermo- 
dynamic potential (TDP) r2(A/, 6, A). In the leading order of the large iVc-expansion it is defined by the following 
expression: 

f d'xn{M,b,A) = -^S,s{a{x),TTa{x)) _ , (11) 

which gives 

\2 , A2 



where 



d^xQiM, 6, A) = J d^xi^L.J!^±^ + _L In [dq][dq] exp(i J d^xqVq^^ , (12) 



qDq = q{Y'^dp + /i7" + VT^f - M)q - A{q^\-lSa) e-^'^'' - A{qdijS^)e^"'^ . (13) 



(Remember, in this formula q is indeed a flavor doublet, i.e. q ~ {qu,qd)'^ ) To proceed, let us introduce in (jl2p - (jl3p 
the new quark doublets, ip Etnd ip, namely: tp — exp{iT3bx)q and ip — gexp(— irafoa;). Since this transformation of 
quark fields does not change the path integral measure in ([T2|) . the expression (fT2|) for the thermodynamic potential 
is easily transformed to the following one: 

d'xn{M, b,A)^J rf2^ (^^~"^o)^ + ^' ^Wj^'il ^'^'^'^^'^^^ '^''K' J '^'^^^^)) ' (14) 



Here and in what follows we will use a rather conventional notation "global" minimum in the sense that among all our numerically 
found local minima the thermodynamical potential takes in their case the lowest value. This does not exclude the possibility that there 
exist other inhomogeneous condensates, different from mOII . which lead to ground states with even lower values of the TDP. 
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where instead of the x— dependent Dirac operator 2? a new x— independent operator D appears 

D = Yid^ - M + ^7" + T3j^b + vT-if - iAri7^ (15) 
The expression (jl4p for the thermodynamic potential is easily transformed to the following one: 

l^(M,5,A)^(^^--°/ + ^%i-^'--^-^-^ 



4G A^c / <P- 



X 



<*'-"'«)^ + '^^iTV,/^,.,I5(rt, ,16, 



4G ' J (27r) 

where the Tr-operation Tisfx stands for the trace in spinor- (s), flavor- (/) as well as two-dimensional coordinate- [x) 
spaces, respectively, and Tr^^ is the respective trace without x— space. Moreover, D{p) =f + ^7° + T-^^^h + vtj,^^ — 
M — i7^ATi is the momentum space representation of the Dirac operator D (jlSp . Obviously, D(p) is a 4x4 matrix in 
the direct product of the spinor and flavor spaces. Since Trs/lnD(p) = IndetD(p), one can evaluate the expression 
(fT6|) with a help of any program of analytical calculations and find 

17(M, 5, A) ^ »""(M, h, A)= ^^' + + i / -^lndet;D(p), (17) 



AG J (27r)2 

where 

det:D(p) = A"* + 2A2(Af2 +p2 ^ ^,2 _ ^2 _ ^2^ 

+ (Af2 + (pi - &)2 -(rj + t.)2) (m2 + (p, + bf - (7? - vf) (18) 

and r] = po + p.. (In order to emphasize the fact that the expression (1171) is divergent, i.e. unrenormalized, we use 
in this TDP notation the superscript "un".) Obviously, the function J7""(M, 6, A) is symmetric with respect to the 
transformation A —A. (At mo = it is also invariant with respect to the M — )• ~M transformation.) Moreover, 
it is invariant under each of the transformations b — >■ —6, /z — >■ — /z and v — > —ly. |f| Hence, without loss of generality, 
we restrict ourselves to the constraints A > 0, > 0, > 0, and v > Q. In the following, we will investigate the 
global minimum point of the TDP ([T7)) just on this region. However, first of all let us consider the case of spatially 
homogeneous condensates, i.e. the 6 = case. 



III. THE CASE OF HOMOGENEOUS CHARGED PION CONDENSATE, 6 = 



Supposing that 6 = in ((T7)) . we obtain after some technical calculations the TDP for the case of spatially 
homogeneous charged pion condensate. 



{M - mo)2 + A^ 
4G 



(Pp 
(2^ 



In- 



{ [{po + pf [Eif] [(po + - [E-^f] }, 



where 



/i//2, E : 



AP. 



(19) 



(20) 



The argument of the ln(a;)-function in (I19|) is proportional to the inverse quark propagator in the energy-momentum 
space representation. Hence, its zeros are the poles of the quark propagator. So, using (jl9p one can find the dispersion 
laws for quasiparticles, i.e. the momentum dependence of the quark {pou, Pad) and antiquark (poa, Pnd) energies, in a 
medium (the full expression of the quark propagator matrix is presented in Appendix B of paper ^2^] ) : 



Pou^E^-fi, pod^E^-p, pou ^ -{EX+ fJ-), Pad^ -{Ea+ fJ-)- 



(21) 



Integrating in over po (see in [2^ for similar integrals), one obtains for the unrenormalized TDP of the system 
at zero temperature the following expression: 



rj""(M, A) 



A2 



4G 



dp 



-{ 



Ea + E. 



+ EX)9{p - eX) + (/i - E^)e{p. ~E-^)], 



(22) 



* Indeed, if simultaneously with b — > —b 01 fi —fi transformations we perform in the integral I I17II the pi — >■ — pi or po — > — po change 
of variables, respectively, then one can easily see that the expression II17I I remains intact. Finally, ii u —i/, we should transform 
pi — > — pi in the integral I I17I1 in order to be convinced that the TDP remains unchanged. 
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where 9{x) is the Heaviside theta-function. It is clear that the TDP (|22p is an ultraviolet divergent quantity, so in 
order to get any physical information one should renormalize it, using a special dependence of such quantities as the 
bare coupling constant G and the bare quark mass mg on the cutoff parameter A (A restricts the integration region 
in the divergent integral in (P^. \pi\ < A). The renormalization procedure for the simplest massive GN model was 
already discussed in the literature, see, e.g., in [3 H^- In a similar way, it is easy to see that, cutting of the 
divergent integral in (|22|) and then using the substitution G = G(A) and toq = mG(A), where 

1 2 /2A\ 

zIMtt (23) 



2G(A) TT V 

and TO, Mq are new free finite renormalization group invariant massive parameter^ (which do not depend on the 
cutoff A), it is possible to obtain in the limit A — >■ c» a finite renormalization group invariant expression for the TDP 
(for details see, e.g., the papers [1^). Namely, 

f7-"(M, A)= VoiM, A) - ^ - ^{e+ + E^- 2^pI+AP + A^ 

+(m - E+)e{^i - E+) + (/i - E^)e{^i - E^)}, (24) 



where 

+ A^ r / )\^2 



Fo(M,A) 



2tt 



(25) 



is a finite renormalization group invariant expression for the TDP (j24p in vacuum, i.e. at /i = and /i/ = 0, taken in 
the chiral limit, i.e. at m = 0. 

In the following, when studying the phase structure, the quantity Mq is still treated as a free parameter, however 
instead of the massive parameter to of the model we will use a dimensionless parameter, a = ttto/Mq. As a result, 
one can see that in the massive NJL2 model the dimensional transmutation effect is absent formally. Indeed, both 
before and after renormalization this massive model is parameterized by one massive- and one dimensionless quantity. 
(Before renormalization the model is characterized by a bare mass toq and a dimensionless bare coupling constant G, 
while after renormalization the mass Mq and dimensionless quantity a are free model parameters.) In contrast, in the 
massless GN-type models, i.e. at toq — 0, the coupling constant G is replaced after renormalization by the massive 
parameter Mq (it is the so-called dimensional transmutation phenomenon). 

In our subsequent calculations throughout the paper the quantity a is fixed by a = ao ~ 0.17. In this case we have 
in the initial NJL2 model the same relation between the pion mass and the dynamical quark mass in vacuum as in 
some NJL-type models in the realistic case of the (3-|-l)-spacetime [26]. 

Investigating the behavior of the global minimum point (whose coordinates are just the gaps M and A) of the TDP 
(l24l) vs chemical potentials, it is possible to establish the phase structure presented in Fig. 1. There, in the phases 
1, 2 and 3 the gap A is vanishing, i.e. these are the normal quark matter phases with a nonzero gap M. However, 
at the boundaries between phases the gap M changes its value by a jump (the details of the investigation, including 
the behavior of gaps, particle densities, meson masses etc, are presented in [2^). For each point {v, fj.) of the vacuum 
region of Fig. 1 we have A = and M « 1.04Mo (the physical meaning of the parameter Mq is described in footnote 
[S]). Finally, one can see in Fig. 1 the homogeneous charged pion condensation phase in which both gaps are not equal 
to zero. What is more interesting for us is that all over this phase the quark number density Uq = —dil.^^"{M, A)/9/i 
is equal to zero, Uq = 0. 

Therefore, in dense (i.e. with nonzero nq) quark matter, mimicked by the initial GN-type model, the phase with 
spatially homogeneous charged pion condensation can not be realized. 



IV. INHOMOGENEOUS ANSATZ FOR THE CHARGED PION CONDENSATE, 6/0 

In this Section, the possibility for the ground state of the initial NJL2 model (1) with spatially inhomogeneous 
charged pion condensate is investigated. We start with the unrenormalized TDP ([T7| which can be rewritten in the 
form 

f]™(M, b, A) = (^^-"^0)^ + ^' + i I _g_ In (^4 ^ ^^2 + Brj + C), (26) 



^ Note, the quantity m is not equal to the physical, or dynamical, quark mass M. The last one is defined by the pole position of the 
quark propagator. Alternatively, it can be found as a gap, i.e. one of the coordinates of the global minimum point of the thermodynamic 
potential. However, parameter Mq is equal to a dynamically generated quark mass in the vacuum and at mo = (a more detailed 
discussion on the physical essence of these parameters is given in l26l ) . 
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FIG. 1. The (/i, v) phase portrait of the model considered 
at T = and v = fii /2 > in the case of spatially ho- 
mogeneous condensates. Numbers 1, 2 and 3 denote dif- 
ferent normal quark matter phases with A = 0, 7^ 0. 
On all the lines of the Figure, first order phase transi- 
tions occur except for the boundary between vacuum and 
homogeneous pion condensation phase, where a second 
order phase transition takes place. Uq is the quark num- 
ber density. 



FIG. 2. The quasiparticle energies (I29p vs pi at /i = 
O.8M0, ly = I.2M0, A = O.35M0, M = O.O6M0, and b = 
0.76Mo. 



where the notation rj ^ po + SiS well as the identity det D(p) = ry* + Ai]^ + Brj + C with 
A = -2(Af2 + b^+pl + iy^ + A^), B = -Spibi^, 

C = (A/2 + b'+pl + u^ + A2)2 - A{p\v^ + b'v^ + + M^v^ + plh'') (27) 
are used. The argument of the In-function in (j26p can be expanded into a product of four hnear multipliers, 

J7""(M, h, A) = ^^'^'"1^^' + ^' + ^ / (S^ [('^ - - "^^^^"^ ~ - % )] , (28) 

where tj^ are presented in Appendix |X] by the expressions (jAl[) . If the quantities A,B and C in ([M]) are defined 
by the relations (j27p . then numerical analysis shows that all the roots rf^ are real valued quantities vs gaps M, A, 6, 
chemical potentials and spatial momentum pi. Moreover, two of the roots, 77^, are negative valued quantities. 

Taking into account the remark after formula (|20p , it is possible to obtain immediately from (|28p the quark-antiquark 
dispersion laws, 

Pou = - M, Pod = vt - Pan =-n~ - ^1, = 77+ - ^. (29) 

Note that at 5 = the quasiparticle energies coincide with the corresponding expressions from (PT|) . At the 
particular values of the chemical potentials and gaps the plots of the quasiparticle energies powv (EH) vs pi are 
presented in Fig. 2. 

Now, it is possible to perform the po-integration in (j28p using the general formula (see Appendix IB)) 

/oo 
dpo hi (po - a) = i7r|a|, (30) 
-00 

where a is a real quantity. As a result, we have 

17«"(M, 6, A) = + - /I ^ [1^ - ''^ I + Ia* - ryri + Im - %+| + Im - % I] . (31) 

To renormalize the expression (1311) we must first regularize it. In this connection, it is necessary to make the following 
remark. In the case of homogeneous condensates (see the previous Section) usually the momentum cutoff regularization 
scheme is used. However it does not work in the case of spatially inhoinogeneous condensates since three-momentum 
is no longer conserve As discussed in the recent papers [13, El, [13, ISl , an adequate regularization scheme in the 



If the momentum cutoff regularization is used in the inhomogeneous case, the TDP acquires some non-physical properties such as 
unboundedness from below with respect to b, etc. As a result, an additional modification of the TDP is needed (for details see in 

[Ulia,!!!]). 
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FIG. 3. Phase portrait of the model when inhomogeneity FIG. 4. Gaps M (line 1), A (Une 2) and inhomogeneity 

of pion condensates is taken into account. Here IPCi,2 wave vector b (line 3) vs fi at fixed f — 1.2A/o. 

denote inhomogeneous charged pion condensation phases, 
Uq is a quark number density, 2 is a normal quark matter 
phase with A = 0, M 7^ (see in Fig. 1). 



case of spatially inhomogeneous phases is that with an energy constraint equal for all quasiparticles. So, dealing with 
spatial inhomogeneity, one can use, e.g., the Schwinger proper-time regularization, dimensional regularization etc. In 
particular, in our recent paper ^35| the symmetric energy cutoff regularization scheme was proposed in considering the 
behavior of chiral density waves in the presence of an external magnetic field in the framework of the four-dimensional 
Nambu-Jona-Lasinio model. There, for each quasiparticle the same (finite) interval of their energy values was allowed 
to contribute to the regularized thermodynamic potential. In the present investigation we will also use the symmetric 
energy cutoff regularization for the TDP, i.e. 

ri"" (M, b, A) — > n'^'-siM, b, A) = (^'^ - ^0? + A' 

- ^ [\Pou\OiA - bonl) + |Pod|0(A - \pod\) + \PoMA - bo^l) + \Podm ~ \Pod\)] : (32) 

where the notations for quasiparticle energies po« etc. are used. Now, let us consider the identity 

f7-^(M, b, A) = (n'-^^iM, b, A) ~ n^^^iM, b, A)l^,^^^^^^^,) + n^^^iM, b, A)l^,^^^^^^^,. (33) 
Clearly, at6 = 0,/i = 0,i^ = for all quantities 775*^2 one finds the relation 17/^*^21 = v^pf~KM^~i-"A2, so that 

^-'(M, b, A)U^__,^^__, ^ ^"""^f - £ ^V^? + M^ + A^ (a - V^f + M^ + A^) . (34) 

Since the expression in parenthesis in (j33p is an ultraviolet (UV) convergent one, i.e. it is a finite quantity in the 
K ^ 00 limit, we see that in (|33|) all the UV divergences are located in the last term which is nothing but energy 
cutoff regularized vacuum thermodynamic potential of the system p4p . Hence, in order to renormalize the TDP 
f2""(Af, b, A) ([3T|) it is sufficient to remove UV divergences from the quantity ([34]) by substituting in (l34l) G = G{A) 
and mo — mG{A) by quantities with an appropriate behavior of G{A) vs A. As a result, we have 

\pou\OiA - \pau\) + \pad\0{A - \pad\) 

+ \po,^0iA - \pou\) + \p,MA - \Pod\) - ^^pI + AP+A^ 6 [a - ^p\ + A/2 + A^ 

where Vo(M, A) is given in (pSj) . 

We have studied numerically the TDP (I55|) as a function of Af , A and b for some physically motivated value of the 
massive parameter m = MqolqI'k^ where ao ~ 0.17 and Mq is a free parameter of the model (see the corresponding 
explanation at the end of the previous section III) . The properties of its global minimum point vs chemical potentials 
give us the phase structure of the model which is presented in Fig. 3. It is easy to see that the phase structure in 
the case of spatially inhomogeneous condensates (see Fig. 3) is the same as in the case of homogeneous ones (see Fig. 



r!""(Af, 6, A) r!-"(A//, 6, A) = Fo(Af, A) ~ ^ - lim \ T ^ 
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FIG. 5. Gaps M (line 1) and A (line 2) vs v at fixed 
= O.SMq. In this case fe = 0. 



FIG. 6. Gaps M (line 1), A (line 2) and inhomogeneity 
wave vector b (line 3) vs v at fixed /i = 0.85Afo. 



1) but with two essential exceptions. Namely, the normal quark matter phases 1 and 3 of Fig. 1 are replaced by two 
inhomogeneous pion condensation phases IPC1.2 in Fig. 3. The behavior of the gaps M and A as well as of the wave 
vector h vs chemical potentials in these phases are presented in Figs 4-6. The TDP ([55)) provides also the expressions 
for the quark number density Uq and isospin density n/, 



917''™(M,6, A) 5r2''^"(M,6, A) 

ni = — 



dm 



(36) 



In particular, as it follows from (1551) and ([55]) 



dp 
Yt\ 
dpi 
2n 



[e{ti - ry+) + 0(m - vi) + e{p^ - vt) + e{pi - %) 
[e{^Ji-lJ+) + e{^i^l^-)\ , 



2] 



(37) 



where the last equality appeares due to the fact that 773 < 0. Using these expressions one can easily prove that in the 
inhomogeneous pion condensation phases IPCi_2 both isospin density nj and quark number density Uq are nonzero 
(see Fig. 3). In contrast, in the homogeneous pion condensation phase of Fig. 3 the density Ug is zero. The behavior 
of ni and rig at some particular values of chemical potentials are presented in Figs 7, 8. Hence, we have proved 
that in the framework of the initial model the charged PC phenomenon with nonzero n^-density is possible only with 
spatially inhomogeneous charged pion condensate. 



SUMMARY AND CONCLUSIONS 



This paper is devoted to the investigation of the so-called charged pion condensation (PC) phenomenon which 
might be observed in dense baryonic matter with different contents of u and d quarks. To simplify the considera- 
tion, we have restricted ourselves to the (l-l~l)-dimensional NJL-type model (1) with quark number fi and isospin /ij 
chemical potentials at zero temperature. Special attention is paid to the influence of spatial inhomogeneity of dif- 
ferent condensates on charged PC phenomenon. Our consideration is performed in the leading order of the large- iVc 
expansion. 

Recall, the charged PC phenomenon was studied recently in the framework of some QCD-like effective theories such 
as NJL models or chiral effective theories in the usual (3-|-l)-dimensional Minkowski spacetime [lol - [l^ . However, the 
existence of the charged PC phase with nonzero baryon or quark number density, denoted below as PCd phase, was 
there predicted without sufBcient certainty. Indeed, for some values of model parameters (the coupling constant G, 
cutoff parameter A etc.) the PCd phase is allowed by NJL models. However, for other physically interesting values of 
G and A the PCd phase is forbidden in the framework of NJL models 11]. Moreover, if the electric charge neutrality 
constraint is imposed, the charged pion condensation phenomenon depends strongly on the bare (current) quark mass 
values. In particular, it turns out that the PCd phase is forbidden in the framework of NJL models if bare quark 
masses reach the physically acceptable values of 5 10 MeV (see in (T^V 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

v/M„ 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 



FIG. 7. Quark number, n,, and isospin, n/, densities vs 
v at )j. = 0.85Mo. 



FIG. 8. Quark number, Uq, and isospin, nj, densities vs 
^ at u = Mq. 



As for investigations of the charged pion condensation phenomenon in the framework of the (l+l)-dimensional 
massive/massless NJL model (1), the results of our recent papers show that the PCd phase is also absent there, if PC 
condensate is spatially homogeneous [H, H^l (see also Sec. Ill of the present paper). However, earlier we have found 
one factor which promotes the creation of the PCd phase at least in (l+l)-dimensional spacetime. It is the finiteness 
of the volume of a physical system |39j . Since such a constraint with certain boundary conditions imposed for any 
system is equivalent to its consideration in a space with nontrivial topology, we have studied in [39l | the initial model 
(1) in the spacetime x (spatial coordinate is compactified) and proved that at some boundary conditions for 
spinor fields the charged PCd phase is realized in the system. 

In this paper, we have proved that a spatial inhomogeneity of PC condensate is also a factor which promotes the 
appearance of PCd phases on the phase diagram of the NJL2 model (1). Indeed, if consideration of dense quark 
matter is performed in terms of homogeneous pion condensates, then PCd phase is absent on the phase diagram of 
the model (1) (see Fig. 1). However, if the spatial modulation of pion condensates is taken into account in the form 
(10), then two PCd phases appear on the phase diagram of the model (1) (those are IPCi and IPC2 phases of Fig. 



In summary, we conclude that charged pion condensation phenomenon of dense and isotopically asymmetric 
quark/hadron matter is more preferable to be spatially inhomogeneous than homogeneous. 

Finally, we would like to discuss the reliability of the main result of our paper and try to predict what might 
happen with the PCd phase in the framework of a more general ansatz for condensates. To simplify the problem, 
let us consider the case of massless NJL2 model (1) with mp = 0, with an evident generalization of the ansatz (|10p . 
Indeed, if toq = 0, then it is possible to study the phase structure of the model in terms of the following simultaneous 
spatial modulations of the chiral and pion condensates 



Evidently, there are three particular cases of i) The choice a = & = correspongs to spatially homogeneous 

condensates and the phase structure of the model for this parameterization was studied in [26] . ii) Then, the choice 
a = is really the ansatz ([T0| at mo = 0. iii) Finally, the phase structure of the model under the constraint & = 
was studied in [30|. Recall that in cases i) and iii) the PCd phase does not appear at the phase diagram. We have 
made preliminary estimations of the phase structure of the massless model (1) in the framework of the ansatz ([55)) 
and found that 1) the phase with M ^ 0, a ^ /S. ^ and & ^ is absent. 2) There is an absolute minimum of 
the TDP corresponding to spatially inhomogeneous PCd phase with A 7^ 0, & 7^ 0, M = 0, a = 0. 3) For the same 
values of chemical potentials there is an equivalent TDP extremum, corresponding to a chiral spiral phase, where 
M 7^ 0, a 7^ 0, A = 0, and 6 = 0. It means that inside an inhomogeneous PCd phase, bubbles of the inhomogeneous 
phase with chiral spiral are allowed to exist and vice versa. By analogy, one might expect that in the more physically 
interesting case with mo 7^ 0, the spatially inhomogeneous charged PCd phase would continue to be present at the 
phase diagram of the model (1), even if an arbitrary more general parameterization of condensates is used. 



3). 



{(t{x)) = M cos(2aa;), (7r3(a;)) = Msin(2ax), (7ri(x)) = Acos(26a;), (7r2(a;)) = Asin(26a;). 



(38) 
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Appendix A: Roots of the equation 77* + Arj^ + B-q + C = 
Using any program of analytical calculations, four roots of this equation can be presented in the following form: 



where 



2A ^/2 R Q 



+ 12c, S ^2A^ + 27B^ - 72AC. (A2) 

Appendix B: Derivation of formula (|30p 

Let us denote the integral in the left hand side of ((30|) by / (recall, there a is a real quantity). 

It well-known that in quantum field theory any loop po -integration is performed in the supposition that po is a 
shorthand notation for po -\- ie ■ sign(po), where e — 0+. In this case the causality of the theory is preserved. Taking 
this circumstance in mind, we see that in (j30l) the integration contour at a > (a < 0) lies above (below) the 
singularity point a — ie- sign(a) of an integrand function. Hence, it is possible to perform in (|30p the Wick rotation of 
the integration contour and to direct it along the imaginary axis of the complex pQ-plane. In thus obtained integral 
one can change an integration variable, po — ^ ipo. As a result, we come to the relation 

/oo /"OO 
dpo In (ipo - a) = i / dpo In (ipo - a) + i / dpn In (ipo - a). (Bl) 
-00 -/ •/ — 00 

In the last integral of fjBip one can again change an integration variable, po — )■ —po. Hence, 

poo pOO nOO 

I='^l dpo\n{ipo - a) + i dpo\n{- ipa - a) = i dpoln [pi + a"^). (B2) 
Jo Jo Jo 



The last integral in (IB2p can be easily taken using the integration by part method. Thus, up to an omitted infinite 
term which does not depend on a, we obtain / = i7r|a|. 
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